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Hunting long-lived gluinos at the Pierre Auger Observatory 
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Eventual signals of split sypersymmetry in cosmic ray physics are analyzed in detail. The study fo- 
cusses particularly on quasi-stable colorless i?-hadrons originating through confinement of long-lived 
gluinos (with quarks, anti-quarks, and gluons) produced in pp collisions at astrophysical sources. 
Because of parton density requirements, the gluino has a momentum which is considerable smaller 
than the energy of the primary proton, and so production of heavy (mass ~ 500 GeV) _R-hadrons 
requires powerful cosmic ray engines able to accelerate particles up to extreme energies, somewhat 
above 10^"^'^ GeV. Using a realistic Monte Carlo simulation with the AIRES engine, we study the 
main characteristics of the air showers triggered when one of these exotic hadrons impinges on a 
stationary nucleon of the Earth atmosphere. We show that i?-hadron air showers present clear 
differences with respect to those initiated by standard particles. We use this shower characteristics 
to construct observables which may be used to distinguish long-lived gluinos at the Pierre Auger 
Observatory. 



I. INTRODUCTION 

There exists "lore" that convinces us that physics be- 
yond Standard Model (SM) should be guided from the 
stabilization of mass hierarchy. The most ubiquitous ex- 
ample is the minimal low energy effective supersymmetric 
theory (MSSM) [1], which requires a scale of supersym- 
metry (SUSY) breaking AgusY ~ 1 TeV to avoid the 
fine tuning problem (/loic ^ ^h/^susy) ^'^^^ the Higgs 
mass {Mh ^ 100 GeV). However, this "naturalness" is 
not favored by prescision tests at colliders, which are con- 
sistent with SM to a great accuracy [2]. Consequently, 
any new physics which may turn on beyond the elec- 
troweak scale needs to be fine tuned at the percent level. 
Moreover, the presence [3] of a tiny, but non- vanishing, 
cosmological constant presents us with a fine tuning prob- 
lem much more severe than the gauge hierarchy problem. 

The solution of last resort to address the cosmological 
constant problem is Weinberg's anthropic approach [4], 
in which there exists an enormous "landscape" of vacua, 
only a small fraction of which have a vacuum energy small 
enough to allow for a natural habitat for observers such 
as ourselves. This approach has been recently rekindled 
by investigations in String Theory which have applied a 
statistical analysis to the large number N of vacua in the 
theory [5] . Among this vast number of metastable vacua, 
there can be a small subset 0(10'*°) exhibiting low scale 
SUSY breaking. Of course the fine tuning required to 
achieve a small cosmological constant implies the need 
of a huge number of vacua, far more than the 0(10'*'') 
characterizing low-scale SUSY breaking. However, the 
density of vacua increases (x Agygy [6]. Therefore, as- 
signing a priori equal probability to each vacuum one 
arrives at a new meassure of fine tuning, which takes 



into account the "entropy" associated with the density 
of vacua, /new ~ -^ff ^susY- Contrary to /lorc require- 
ments, /new clearly favors a large SUSY breaking scale. 
For example, for AgusY ~ 10^° GeV, ©(lO^o") vacua 
become available, enough to fine tune both the cosmo- 
logical constant and the Higgs mass. If we live in this 
neighborhood of the landscape, collider data would be 
expected to point to the SM rather than SUSY. How- 
ever, one pays a price for throwing away MSSM, since it 
provides a potential explanation for both dark matter [7] 
and the LEP results favoring the unification of the three 
SM gauge couplings [8]. 

Split SUSY [9] is a relatively new variant of SUSY 
which may facilitate the required fine tuning and simul- 
taneously preserves the achievements of the MSSM. In 
this model the bosonic superpartners are heavy, while the 
extra fermions retain TeV-scale masses thanks to protec- 
tion by chiral symmetry. Although split SUSY does not 
provide a dynamical explanation for the hierarchy prob- 
lem, the assumption of a large-scale SUSY breaking leads 
to important information on the underlying parameters 
and on measurable physical quantities [10]. In particular, 
analyses of one loop [11] and two loops [12] running of 
the RG equations, show that split SUSY preserves uni- 
fication of couplings. Additionally, as in the MSSM, the 
lightest supersymmetric particle provides a possible can- 
didate for cold dark matter [13]. 

It is clear that split SUSY opens new territory for 
model builders, gauginos have a symmetry that protects 
their masses, namely the i?-symmetry, so building models 
where scalars are very massive is quite natural in theo- 
ries where this symmetry is not broken, for example in D- 
term breaking models; it can also happen in theories with 
extended supersymmetries and there are already several 
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papers with string inspired models of split SUSY [14]. 
Hence there is a strong motivation for phenomenologi- 
cal studies, including implications for collider-based mea- 
surements [15-17], electric dipole moments [18], Higgs 
physics and electroweak symmetry breaking [19], and cos- 
mic ray physics [16, 20, 21]. The latter is the main focus 
of the present study. 

An intriguing prediction of split SUSY, which repre- 
sents a radical departure from the MSSM, is the longevity 
of the gluino. As mentioned above, in split SUSY the 
squarks are very massive and so gluino decay via virtual 
squarks becomes strongly suppressed, yielding a g life- 
time of the order of [9] 



(1) 



where Mg is the gluino mass. Quasi-stable colorless R- 
hadrons (i.e., carrying one unit of i?-parity) are expected 
to be born when such long-live gluinos become confined 
with quarks, anti-quarks, and gluons [22]. 

Very strong limits on heavy isotope abundance in turn 
require the gluino to decay on Gyr time scales [23] , lead- 
ing to an upper bound for the scale of SUSY break- 
ing 0(10^'^) GeV. More restrictive bounds on Asusy 
can be determined from cosmological considerations [24] . 
Specifically, gluino decays would disturb predictions of 
big bang nucleosynthesis (BBN), or distort the cos- 
mic microwave background (CMB), or alter the diffuse 
gamma ray background. The details depends on both, 
Mg and Tg. For example, for 1 < Mg/TeV < 5, to avoid 
altering the abundances of D and ^Li, Tg < 100 s, im- 
plying Asusy < 10^" GeV [24]. The relic abundance 
of lighter gluinos, Mg < 500 GeV, is constrained by 
COBE [25], WMAP [26] and EGRET observations [27]. 
On the one hand, gluinos that decay during or after the 
thermalization epoch can distort the CMB spectrum [28] 
and so are limited by COBE/WMAP observations. On 
the other hand, gluinos that decay after the recombina- 
tion epoch give rise to pions which subsequently decay 
into 7-rays that free-stream to us. The contribution of 
such a decay chain to the diffuse 7-ray background [29] 
is limited by EGRET observations. 

Long-lived gluinos are also constrained by collider 
searches. Charged i?-hadrons can be observed as they 
cross the detector either by their time delay relative to 
ultrarelativistic particles [30], or by their anomalously 
high ionization energy loss [31]. Besides the energy de- 
position of neutral i?-hadrons in the calorimeter is rather 
soft, and so when they are produced in association with 
a high-pT jet they can be observed in the monojet chan- 
nel + missing energy ■ CDF Run I data [32] found a 
bound of Mg > 170 GeV [16]. In addition, i?-hadrons 
can become stopped gluinos by losing all of their mo- 
mentum and coming to rest in the calorimeter [17]. The 
D0 Collaboration [33] has recently searched for stopped 
gluinos decaying into a single jet and a neutralino. The 
non-observation of monojets (above the expected back- 
ground from cosmic-muon induced showers) in Run II 
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FIG. 1: Limits on long-lived gluinos. The cross-hatched bands 
indicate excluded regions of the Mg-Tg plane from anomalous 
heavy isotope abundance [23], CMB [25] and EGRET [27] ob- 
servations, BBN predictions [24], and collider data [16, 33]. 
Contours of constant values of Asusy are also shown by 
solid (lO^" GeV), dash-dotted (lO^^ GeV), and dashed lines 
(10^^ GeV). 



data implies Mg > 270 GeV for rg < 3 hr. All these lim- 
its are shown in Fig. 1. As we will show here, the study 
of hadronizcd gluinos originating in distant astrophysical 
sources, provide a viable experimental handle in the re- 
gion 300 < A/g/GeV < 500 - 10^ < Tg/yr < 10^ which 
is yet unexplored. 

The main goal of this paper is to describe a full-blown 
Monte Carlo simulation of R air showers, and uncover 
observables which may be exploited by new experiments 
like the Pierre Auger Observatory [34]. This analysis 
expands on previous work [35] by including all possible 
R interactions and analyzing in detail the potential of 
the surface array. Before describing the simulation, we 
introduce in the following section the main properties of 
i?-hadron interactions 



II. COSMIC R's 

The origin cosmic rays is still an open question, with 
the degree of uncertainty increasing with rising en- 
ergy [36] . Theoretically, one expects the cosmic ray spec- 
trum to fall off somewhat above 10^"'^ GeV, because 
the particle's energy gets degraded through interactions 
with the cosmic microwave (protons and nuclei) and ra- 
dio (photons) backgrounds, a phenomenon known as the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff [37]. The most 
recent data from the Pierre Auger Observatory in fact 



3 



do not indicate (yet) any excess beyond the expected 
cutoff [38]. Because of tlie rapid energy degradation, 
the maximum attainable energy in far away sources can 
be considerably higher than the GZK-threshold. How- 
ever, measurements of the GeV diffuse gamma ray flux 
significantly constrain the cosmic ray production inte- 
grated over redshift, and consequently limit the maxi- 
mum energy of these particles. Specifically, the interme- 
diate state of the reaction P7cmb ^ -^^r is dominated 
by the A+ resonance (because the neutron decay length 
is smaller than the nucleon mean free path on the relic 
photons). Hence, there is roughly an equal number of 
TT^ and tt". Gamma rays, produced via tt*^ decay, subse- 
quently cascade electromagnetically on the cosmic radi- 
ation fields through e+e" production followed by inverse 
Compton scattering. The net result is a pile up of gamma 
rays at GeV energies, just below the threshold for further 
pair production. Therefore, if the distribution of cosmic 
ray sources is homogeneous and each source is character- 
ized by a hard injection spectrum oc E~^, then EGRET 
measurements in the 100 MeV - 100 GeV region [27] limit 
the maximum proton energy ~ 10^^'^ GeV [39]. Since 
Fermi's acceleration mechanism predicts a rather steeper 
spectrum cx E~'^ [36], one can assume a maximum proton 
energy El%,^ ~ IQi^-^ GeV. 

Gluinos are flavor singlets of a color SU{3) octect 
that interact strongly with the octect of gluons and can 
combine with quarks, anti-quarks, and gluons to form 
colorless hadrons [22]. The bosonic states, gqqq, are 
generically called i?-baryons, whereas the fermionic sates, 
gqq and gg, are called i?-mcsons and i?-glueballs, re- 
spectively. Very little is certain about the spectroscopy 
of these strongly interacting particles. The most rel- 
evant feature is (perhaps) the difference in mass be- 
tween i?-mesons (-glueballs) and i?-baryons, because if 
^Rm(g) + nT'N > Mji^ + TO^, then there are exothermic 
conversions of Rmig) into Rb as the i?-hadrons propagate 
in the atmosphere. 

-R-hadrons states should be produced in pairs through 
pp collisions at powerful cosmic ray engines (e.g., pro- 
tons undergoing acceleration in compact jets of relativis- 
tic plasma interact with those in the surromiding gas). 
The average energy of the produced R in the target sys- 
tem is 
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2m„ ^ 



(2) 



where is the energy of the proton undergoing ac- 
celeration at the source, E^j^ ~ yi/2 is the average 
R energy in the center-of-mass (cm.) of the pp col- 
lision, and s is the square of the energy in the cm. 
of the parton-parton collision. Now, by restricting R- 
production to large cm. energies (say, s > 16M^), from 
Eq. (2) we obtain the maximum energy of cosmic R^s, 
E^^^ < lO'^Mfl. Thus, the R spectrum cuts off at lower 
energy than the cosmic ray spectrum. Since these par- 
ticles originate from cosmological distance d, to reach 



the Earth the gluinos must be remarkably long lived 
Tg > 100 (Mg/500 GeV) (d/Gpc) yr. Moreover, to avoid 
deflections on the extra-galactic magnetic fields and the 
consequent energy loss due to pair production and other 
mechanisms (such as synchrotron or bremsstrahlung ra- 
diation), the i?-hadron has to be neutral [40]. The overall 
intensity of i?-hadrons is constrained by its accompany- 
ing pion flux, which decays into 7-rays and neutrinos that 
can be confronted with existing data [27, 41]. Unfor- 
tunately, the expected flux of ultra-relativistic (Lorentz 
factors ~ 10^) ii-hadrons is found to be very low (less 
than 6 particles per km^ per millenium [20]), and so the 
only experimental method with potential is observation 
of their interactions in large volumes of the Earth's at- 
mosphere. 

When a hadronized gluino impinges on a stationary 
nucleon of the Earth atmosphere a large number (over 
140 when summed over all i?-hadrons) of scattering pro- 
ceses are possible [42]. Interactions of i?- meson states 
include: (i) 2 ^ 2 processes, such as purely elastic, 
(e.g. gdd + uud gdd + uud), charge exchange (e.g. 
gdd + uud gud + udd), and baryon exchange (e.g. 
gdd + uud gudd + ud); (ii) 2 X processes in- 
cluding normal inelastic scattering (e.g. gdd + uud 
gud + udd + dd) and inelastic scattering with baryon ex- 
change (e.g. gdd + uud guud + ud + du + dd). Since 
the flnal-state pion is so light, processes with baryon ex- 
change would be kinematically favored. However, these 
processes could be dynamically suppressed because the 
exchange of two quarks is required. Interaction of R- 
baryon states include purely elastic, charged exchange, 
and normal inelastic scattering. No baryon exchange is 
possible because of the negligible probability for a gqqq 
to interact with a pion in the nucleus. Furthermore, this 
process would be kinematically strongly disfavored. Con- 
sequently, i?-mesons can convert into i?-baryons, but not 
vice versa. Interactions of i?-glueballs are expected to be 
similar to those of i?-mesons. This is because a, g is able 
to split into a qq state, suggesting that a gg interacts like 
(and mix(« with) gqq states. 

To establish which of these processes dominates, aside 
from a model describing the target (neutron or proton), 
the relative couplings of all the processes must be known. 
The latter requires the calculation of the Clebsh-Gordon 
coefficients of isospin-related processes, and the evalua- 
tion of all adittional dynamical effects for all processes. 
To parametrize our ignorance about QCD interactions, 
here we will consider all the relevant processes mentioned 
above (5 for i?- mesons, and 3 for i?-baryons), assigning 
them different probabilities ad hoc so as to explore the 
entire parameter space. 

Predicting the total cross section of an ii-hadron scat- 
tering off a nucleon is non-trivial. However, because of 
the high cm. energies imder consideration in this pa- 
per, the cross section can be safely approximated by the 
geometrical cross section. Moreover, since the size of 
the i?-hadron is roughly the same as the size of the ac- 
companying hadron system, the total cross section for 
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nucleoli scattering can be approximated by the asymp- 
totic values for the cross sections for normal liadron scat- 
tering off nucleons. Therefore, for i?-baryons we take 
aR,-p{^ ~ 10^ GeV) « 140 mb [43]. This corre- 
sponds to a cross section for scattering off air molecules 
(^Rb-airi^/s — 10^) GcV » 520 mb, yielding a mean free 
path in the atmosphere Xr,^ = Wair/f/j^-air ~ 47 g/cm^, 
where we have taken rriair ~ 2.43 x 10~^^ g (correspond- 
ing to an atomic mixture of 78% N, 21.05% O, 0.47% 
Ar and 0.03% of other elements). At this cm. en- 
ergy, the TT-air cross section is roughly 90% of the p- 
air cross section [44], hence for i?-mcson states we set 
Xr^ = 52 g/cm^. i?-glueballs are expected to have the 
same cross section as i?-mesons. This is because the ge- 
ometrical cross section is approximated by the high en- 
ergy hadron cross section, where gluon exchange domi- 
nates (the gg coupling is a factor 9/4 larger than the qg 
coupling, but a meson has 2 quarks, resulting in a cross 
section of a gg state which is (9/4)/(l -|- 1) ~ 1 times the 
cross section for a gqq state). With this in mind we set 
Aii„ ^ XRg . 

In analogy to a billard ball moving through a sea of 

ping-pong balls, the R suffers very little energy loss as 
it traverses the atmosphere. Then for i?-flippers (i.e., 
Rm(g) Rb) we assume that the emitted pion has an 
energy E-„ » Yrn^Tx, where V is the Lorentz factor of the 
incoming i?-hadron. This means that most of the energy 
~ TMg is carried by the accompanying i?t produced in 
the interaction. Following [20], for inelastic collisions we 
parameterize the fractional energy loss per collision as 
-f^inci ~ (A^k/GcV)"^. For completness, a derivation of 
this relation is given in the Appendix. 



III. THE PIERRE AUGER OBSERVATORY 

The Pierre Auger Observatory (or simply 
"Auger") [34] is designed to study cosmic rays with en- 
ergies above about 10^ GeV, with the aim of uncovering 
their origins and nature. Such events are too rare to 
be directly detected, but the direction, energy, and to 
some extent the chemical composition of the primary 
particle can be inferred from the cascade of secondary 
particles induced when the primary impinges on the 
upper atmosphere. These cascades, or air showers, have 
been studied by measuring the nitrogen fluorescence 
they produce in the atmosphere or by directly sampling 
shower particles at ground level. Auger is a hybrid 
detector, exploiting both of these well established 
techniques, by employing an array of water Cerenkov 
detectors overlooked by fluorescence telescopes. On 
clear moonless nights, air showers are simultaneously 
observed by both types of detectors, facilitating powerful 
reconstruction methods and control of the systematic 
errors which have plagued cosmic ray experiments to 
date. 

The observatory is now operational on an elevated 
plane in Western Argentina and is in the process of grow- 



ing to its final size of 3,000 km^. The surface detector 
(SD) consists of an array of 1600 water tanks deployed on 
an hexagonal grid with spacing of 1.5 km. These tanks 
detect the Cerenkov light produced by shower particles 
crossing their 1.2 mx 10 w? water volume, thanks to three 
9-inch photo-multipliers. The fluorescence detector (FD) 
consists of four ensambles of six telescopes, each of which 
has a field of view of 30° vertically and 30° horizontally 
(i.e., 180° for each fluorescence detector site). The geog- 
raphy of Northern site would accommodate a larger array 
(of up to 10,370 km^ [46]), allowing higher sensitivity to 
the low flux of cosmic i?'s. 

Identifying showers themselves is usually straightfor- 
ward, as there is essentially no "background" for the de- 
tectors, at least above their energy threshold. In the case 
of Auger, the threshold for the surface detector is around 
10^'^ GeV, below which less than 10% of the showers can 
trigger three tanks or more, as required. However, full de- 
tection efficiency (i.e., 100% or "saturated acceptance") 
is achieved only around 10^'^ GeV for showers with zenith 
angle lower than 60°, and lower energy showers are usu- 
ally discarded to avoid any complication caused by the 
energy dependence of both the detection efficiency and 
the energy resolution. For fluorescence detectors, showers 
with energies as low 10^ GeV can be observed. However, 
the corresponding acceptance is rc^latively low, since the 
total intensity of the fluorescence light does not allow de- 
tection from a large distance, and the shower maximum 
is then usually above the field of view of the telescopes, 
which prevents accurate reconstruction. Like for any flu- 
orescence detector, the acceptance of the eyes of Auger 
increases with energy (as bigger showers can be seen from 
larger distances) and depends on the atmospheric condi- 
tions. However, a precise determination of the fluores- 
cence detector acceptance is not crucial, thanks to its hy- 
brid nature, the energy differential flux (or "spectrum") 
is not obtained from the fluorescence detector, but from 
the surface detector whose absolute acceptance is essen- 
tially geometrical above saturation and thus is controlled 
within a few percent at most. 

Identifying the primary particle species is somewhat 
more difficult as one has to search for differences in the 
shower development, which are usually relatively small 
and subject to fluctuations associated with the stochas- 
ticity of the first interactions [45]. However, as we discuss 
in the next sections the showers initiated by i?-hadrons 
have very distinctive characteristics and can be easily iso- 
lated from background. 



IV. AIR SHOWER SIMULATIONS 

The AIRES simulation engine [47] provides full space- 
time particle propagation in a realistic environment, tak- 
ing into account the characteristics of the atmospheric 
density profile (using the standard US atmosphere [48]), 
the Earth's curvature, and the geomagnetic field (calcu- 
lated for the location of Auger with an uncertainty of a 



5 



few percent [49]). 

The following particles are taken into account in 
the AIRES simulations: photons, electrons, positrons, 
muons. pions, kaons. eta mesons, lambda baryons, micle- 
ons, antinucleons, and nuclei up to = 36. Nucleus- 
nucleus, hadron-nucleus, and photon-nucleus inelastic 
collisions with significant cross-sections arc taken into ac- 
count in the simulation. The hadronic processes arc sim- 
ulated using different models, accordingly to the energy: 
high energy collisions are processed invoking an external 
package (sibyll 2.1 [50] or QGSJEtII [52]), while low en- 
ergy ones arc processed using an extension of the Hillas 
splitting algorithm (EHSA) [53]. The threshold energies 
separating the low and high energy regimes used in our 
simulations are 200 GeV and 80 GeV for the sibyll and 
QGSJET, respectively. The EHSA low energy hadronic 
model used in AIRES is a very fast procedure, effec- 
tively emulating the major characteristics of low energy 
hadronic collisions. The model is adjusted to retrieve 
similar results as the high energy hadronic model for 
energies near the transition thresholds previously men- 
tioned, and the low energy cross sections are calculated 
from parameterizations of experimental data. A com- 
plete discussion on the low energy hadronic models is 
clearly beyond the scope of this paper. A separate re- 
port on this subject will be published elsewhere [54]. 

The AIRES program consists of various interacting 
procedures that operate on a data set with a variable 
number of records. Several data arrays (or stacks) are 
defined. Every record within any of these stacks is a par- 
ticle entry and represents a physical particle. The data 
contained in every record are related to the character- 
istics of the corresponding particle. The particles can 
move inside a volume within the atmosphere where the 
shower takes place. This volume is limited by the ground, 
the injection surfaces, and by vertical planes which limit 
the region of interest. Before starting the simulation all 
the stacks are empty. The first action is to add the first 
stack entry, which corresponds to the primary particle. 
Then the stack processing loop begins. The primary is 
initially located at the injection surface, and its down- 
wards direction of motion defines the shower axis. After 
the primary's fate has been decided, the corresponding 
interaction begins to be processed. The latter generally 
involves the creation of new particles which are stored 
in the empty stacks and remain waiting to be processed. 
Particles entries are removed when one of the following 
events happen: (a) the energy of the particle is below 
the selected cut energy; (6) the particle reaches ground 
level; (c) a particle going upwards reaches the injection 
surface; (d) a particle with quasi horizontal motion ex- 
ists the region of interest. After having scanned all the 
stacks, it is checked whether or not there are new particle 
entries pending further processing. If the answer is posi- 
tive, then all the stacks are scanned once more; otherwise 
the simulation of the shower is complete. 

AIRES has been successfully used to study several 
characteristics of high energy showers, including com- 



TABLE I: Interaction probabilities for the possible Ji-hadron 

scattering processes. 



Hadron 


Pi 


P2 
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Pi 


Pi 


Rb 
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0.1 
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0.05 
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parisons between hadronic models [55], influence of the 
LPM effect [56], muon brcmsstrahlung [57], and geomag- 
netic deflections [49] on the shower development. AIRES 
has been also successfully used to determine the effi- 
ciency of Auger for quasi-horizontal showers generated 
by r-neutrinos [58], to estimate the flux of atmospheric 
muons [59] , and to study the production of black holes in 
TeV-scale gravity models [60]. For the present analysis, 
we prepared a new module to account for the simulation 
of cosmic i?'s. The module includes external parameters 
such as the type of primary hadron [Ri,, or Rm), its mass 
Mr, its charge, and its primary energy . We adopt 
the atmospheric mean free path derived in the previous 
section. 

The total interaction probability is managed by five pa- 
rameters (Pj, i = 1, . . . , 5) which measure the weights of 
the different processes. (1) The parameter Pi measures 
the probability of an elastic scattering. The program 
emulates this process by transferring a small part of the 
R energy 1 TcV) to an air nucleus which is injected 
into the shower. (2) The parameter P2 measures the 
probability of an clastic scattering with charge exchange, 
in which there is again a small transfer of energy to an 
air nucleus, but now the R primary also flips its charge. 
(3) The parameter P3 accounts for baryon exchange. In 
this case, a pion with energy = mTr{E^^ /Mr) is in- 
jected into the shower evolution, where E^^ is the en- 
ergy of the P-hadron before the collision. If the in- 
coming P-hadron is neutral, then the outgoing pion is 
charged and vice versa. The P-hadron mutates into a 

baryon with E^^ = E^^ — E^^ ■ (4) In the inelastic pro- 
cess, controlled by P4, the P-hadron transfers an energy 
EcoW = E^^^Kinei to the shower. The process is simulated 
with the help of the standard packages of AIRES for nu- 
cleon and pion collisions. A neutron (of energy Pcou) 
scatters off an air nucleus if the primary is a P^, and a 
proton if it is a P^. If it is a Pj„ the projectile particle 
(simulating the P collision) is a tt", whereas if it is a P^ 
the projectile is a tt^. All the secondaries resulting from 
this interaction are considered in the subsequent shower 
evolution. Finally, (5) P5 controls the inelastic scattering 
with baryon exchange. It is simulated as a combination 
of processes (3) and (4), i.e, the emission of a pion fol- 
lowed by an inelstic collision. In our simulations we take, 
P™ = 0.8 - P4™, with < P4™ < 0.8. 

For the simulation engine, the shower starts when the 
P-hadron is added to the previously empty stack. The 
injection surface is located at the top of the atmosphere, 
spacing the interaction point according to an exponential 
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FIG. 2: Energy fraction deposited in the atmosphere during i?-hadron air showers as a function of the zenith angle Q (left) 
and sec0 (right). The curves represent an average over the different species. We note that the distinguishing power among 
the various species is very limited for the entire range of P™. 



distribution with mean equal to Xr^^^^-^-^ ■ The interaction 
probabilities Pi are given in Table I. The i?-hadron is 
tracked until it reaches ground level or else its energy 
is degraded below 100 MeV. New stack entries are ap- 
pended to the existing lists for every SM particle pro- 
duced in the i?-interactions. These entries are then re- 
peatedly processed sequentially by means of the algo- 
rithms implemented in AIRES. 

As an illustration, we have run a set of air shower simu- 
lations, with Mr = 500 GeV, and E^g^ = lO^-^ GeV [61]. 
In Fig 2 we show the energy fraction dissipated into "vis- 
ible" particles in i?-hadron air showers, as predicted by 
our simulations. One can see from the figure that there 
needs to be sufficient pathlength for the R, with its low 
inelasticity, to lose sufficient energy. The experimen- 
tally interesting region to search for i?-hadrons is then 
70° < 6 < 90°. 

Because of the very low inelasticity of i?-air inter- 
actions the leading particle retains most of its energy 
all the way to the ground, while the secondary parti- 
cles promptly cascade to low energies as for any other 
air shower. This results in an ensemble of mini-showers 
strung along the trajectory of the leading particle. Since 
the typical distance between mini-showers is about 10 
times smaller than the extent of a single longitudinal 
profile, it is not possible to resolve the individual showers 
experimentally. Instead one observes a smooth envelope 
encompassing all the mini-showers, which extends from 
the first interaction all the way to the ground, see Fig. 1 
in Ref. [35]. The i?-hadron air showers then present a 
distinct profile: the flatness of the longitudinal develop- 



TABLE II: Primary proton energy Ej^ required to produce 
the same total signal at ground level than a _R-hadron with 
^lab ^ -^q9.7 Qgy^ different zenith angles. 



e 


60° 


66° 


72° 75.5° 


78° 


81.4° 


84° 


87.1° 


El"-^ (GeV) 


10^" 




^09.3 ^(,9.5 


^(,9.6 


^q9.8 







ment is unique to the extremely low inelasticity of the 
scattering, and can be easily isolated from background. 
However, it turns out that there is a sharp cutoff in the 
production of cosmic i?'s at E^^^ « 10^-^ GeV [20], which 
unfortunately leads to showers below detection threshold 
for the fluorescence method (except for a very small aper- 
ture comprised of regions close to the telescope). 

In Fig. 3 we show a comparison of the total signal 
at ground level for R and proton air showers. To com- 
pute the total signal of a single shower, we flrst consider 
all particles reaching the ground with a distance to the 
shower axis r > tq, with tq conveniently chosen to be 
250 m, and for each particle we simulate the Cerenkov 
detector response. The total signal of a given shower, 
S{r > To), is the sum of each particle's individual sig- 
nal normalized to 1 for proton showers of 10^" GeV and 
incident zenith angle of 60°. It is clear that the total 
R signal at ground level increases with zenith angle, be- 
cause of the larger slant-depth. This is in sharp contrast 
to proton showers, in which the signal is reduced with 
increasing because of the greater shower age. 

As discuss in Sec. Ill, the relation between the signal 
observed at the surface detectors and the primary energy 
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FIG. 3: Average total signal beyond 250 m from the shower 
core versus shower zenith angle. Signals are normalized to 1 
for the case of lO^'' GeV proton showers inclined 60°. The 
dashed lines correspond to proton showers of primary energy: 
(a) 10^ ° GeV, (b) lO^'^ GeV, (c) 10"' ° GeV, (d) W^'> '' GeV, 
(e) 10^^" GeV. The solid line corresponds to i?-hadron show- 
ers. The lines are only to guide the eye. 

is determined using hybrid events in which the fluores- 
cence eyes are thought to provide a reliable measurement 
of the total energy. For proton showers the surface detec- 
tors sample about 1% to 10% of the shower energy. Be- 
cause of the electromagnetic component recycling, the R- 
hadron produces a somewhat larger signal at ground level 
than one would expect from standard baryonic showers. 
As one can check in Table II, for large zenith angles if one 
assumes the shower properties are the characteristics of 
proton showers then the total primary energy would be 
overestimated. Note that this aspect is not compensated 
by the calibration procedure, because the i?-component 
does not deposit significant energy in the region of the at- 
mosphere used in the fluorescence-based calibration. In 
summary, although the total contribution to the shower 
energy is small, the i?-hadron deposits a disproportional 
large fraction of their energy close to the ground. Con- 
sequently, cosmic i?'s would induce a significant signal 
in the surface array but not in the fiuorescence eyes. In 
what follows we use these shower characteristics to con- 
struct observables which may be used to distinguish R- 
hadron from traditional cosmic ray showers. 

V. GROUND ARRAY SIGNAL 

The surface detectors of the Auger Observatory are ca- 
pable of measuring the signal associated to an incoming 




2 4 6 S 10 12 14 

Delay time [pa] 



FIG. 4: Lateral vs. arrival time delay distribution for a ver- 
tical shower initiated by a 10^'^ GeV proton. 



shower as a function of time. Since a high energy event 
triggers many detectors, placed at different distances 
from the shower axis, it is possible to reconstruct the 
lateral-time distribution of the signal S{r,t). S{r,t)dt 
gives the amount of signal at a (3-dimensional) distance 
r from the shower axis, produced at the time interval 
[t, t + dt\. For convenience, the origin of times is defined 
for each point in the ground surface as the instant where 
a plane orthogonal to the shower axis, sinchronized with 
the primary particle and moving towards the ground at 
the speed of light, intersects the corresponding point. In 
this way, S{r, t) is necessarily zero for negative times. 
With this definition, the time t is frequently called "ar- 
rival time delay." 

The total signal at a given distance from the shower 
axis is the signal accumulated over all times, that is, 

/•OC 

Stotir) = / S{r,t)dt . (3) 
Jo 

Other quantities that are usually used in the analysis of 
SD signals are: 

• Shower front arrival time, Iq {Iq > 0). This is the 
time corresponding to the first nonzero shower sig- 
nal at the given point, to is directly related to the 
shower front curvature. 

• Partial rise times, t^, defined as the time elapsed 
until the accumulated signal is a fraction x of the 
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FIG. 5: Same as figure 4 but in the case of 10^'^ GeV protons 
inclined 75 degrees. 

total signal, that is, 




Common values of x are: 10%, 50%, and 90%. tx 
is a growing function of r, especially far from the 
shower axis. 

The lateral-time distribution of the signal is a SD ob- 
servable capable of characterizing showers initiated by 
cosmic rays. Consider, for example a typical shower ini- 
tiated by a vertical proton. In Fig. 4 the corresponding 
lateral-time signal distribution is displayed using a false 
color (or grayscale) diagram. From this figure it is pos- 
sible to notice the main features of such a distribution: 
(i) most of the shower particles arrive near the shower 
axis, that is, the signal lateral distribution Stotir) de- 
creases with r; (ii) tQ{r) increases with r, as expected, be- 
cause particles must travel longer distances and undergo 
more interactions to reach positions located far from the 
shower axis; (iii) the time interval of the signal at a given 
point grows with r. In the example of Fig. 4 it goes from 
some 4 /xs at r = 300 m to about 12 fis for r > 3000 m. 

If the inclination of a shower is increased, the thick- 
ness of the air layer placed between the point where 
the cosmic particle enters the atmosphere and ground 
level, also increases. As a result, the age of the detected 
shower increases too. In the case of showers initiated by 
hadronic primaries like protons and nuclei, the aging of 
inclined showers at ground becomes evident for inclina- 
tions larger than 65 degrees, because of the practically 




Delay time [pa] 



FIG. 6: Same as figure 4 but for showers initiated by 10^'^ 
GeV J?-hadrons inclined 75 degrees. 



complete attenuation of the electromagnetic component 
of the shower. For such inclinations, the muonic com- 
ponent becomes very important (see Fig. 2 of Ref. [49]), 
because it produces significant modifications in the de- 
tected signal. In particular, the shower front becomes 
flatter, and the signal is concentrated within a relatively 
small time span. These caracteristics show up clearly in 
Fig. 5, where S{r,t) is plotted for showers initiated by 
protons with the same energy than the showers shown in 
Fig. 4, but for an inclination of 75 degrees with respect 
to the vertical. 

The lateral-time signal distributions of inclined show- 
ers initiated by i?-hadrons present a substantially dif- 
ferent aspect, when compared with the proton case. In 
Fig. 6 the signal distribution corresponding to 10^'^ GeV 
i?-hadron showers inclined 75 degrees, is displayed. A 
comparison with the distribution of Fig. 5 leads to the 
following conclusions: (i) the i?-hadron distribution is 
slighlty more concentrated near the shower axis, and (ii) 
the time span of the signal is substantially larger than 
in the proton case (note that the primary energy of the 
proton showers has been chosen accordingly with Table II 
such that the amount of signal for r > tq is, on average, 
the same for both primaries). 

The last feature of the lateral-time distribution of R- 
hadron showers is certainly the most clear signature of 
such events that could be found in our simulation study. 
Combined, in the case of hybrid events, with a neatly 
different longitudinal development, and inconsistent en- 
ergy measurements, i?-hadron events can be clearly dis- 
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FIG. 7: tio — to, tso — to, and t^o — to plotted versus the 
logarithm of the distance to the shower axis. The data corre- 
sponds to lO'''^ GeV protons inclined 75 degrees. 



tinguished from hadronic ones, and also from neutrino 
initiated showers where the FD energy determination will 
be very different from the present case of i?-hadron show- 
ers. 

The different time span of signals can be quantified 
more precisely studying the observables iio, iso, and tgo- 
Figures 7 and 8 contain plots of these observables as func- 
tions of log]^Q(r), in the case of proton and i?-hadron 
showers, respectively. The larger time span of the sig- 
nals in the i?-hadron case is evident for all the plotted 
quantities (note the different time scales used in either 
figure) . 

The total signal as a function of the distance to the 
shower axis, Stot(r), called lateral distribution, is other 
fundamental observable that can be measured with the 
Auger surface detectors. It is the most important observ- 
able for SD energy determination in the case of showers 
with an inclination of less than 65 degrees. 

We have studied the behaviour of the lateral distribu- 
tion in the case of i?-hadrons, comparing the results with 
the corresponding ones for showers initiated with pro- 
tons. The results, displayed in Fig. 9, clearly show that 
the signal corresponding to i?-hadron showers is more 
concentrated near the shower axis. This also implies a 
different slope for the distributions that could eventually 
be measured. As a reference, a qualitative indicator of 
the signal threshold of Auger detectors is also shown in 
the figure (dotted line). 
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FIG. 8: Same as figure 7, but for 10®'^ GeV i?-hadrons in- 
clined 75 degrees. 
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FIG. 9: Water Cherenkov signal plotted versus the logarithm 
of the distance to the shower axis. The solid (dashed) his- 
togram correspond to Rhadron (proton) primaries. The hor- 
izontal dotted line indicates (approximately) the threshold of 
Auger surface detectors. 
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VI. CONCLUSIONS 

In this paper we have analyzed the characteristics of 
air showers produced by ghiino-containing hadrons. Us- 
ing AIRES, we have performed a high statistics set of full 
quality showers initiated by i?-hadrons. We have consid- 
ered both i?-baryon and i?-mcson primaries. The analy- 
sis of standard observables that can be measured by hy- 
brid air shower experiments like Auger shows that atmo- 
spheric cascades initiated by i?-hadrons are significantly 
different from "classical" showers, such as for example, 
those initiated by cosmic protons. Our study indicates 
that if cosmic J?-hadrons do exist they would produce 
a particular signature that will be visible at Auger: at 
ground level, the i?,-hadron showers are characterized by 
the presence of a strong electromagnetic component at 
all zenith angles. This implies a much longer time span 
for the signal, in comparison with proton showers. At the 
same time, the minishowers generated by the passage of 
the i?-hadron across the atmosphere produce narrower 
lateral distributions than the corresponding ones for the 
proton case. 

If i?-hadron events are analyzed with the; standard pro- 
tocol for hadronic primary showers, a series of inconsis- 
tencies will be present. In particular, the energy de- 
termination via ground signal analysis of very inclined 
showers [62] will likely lead to a primary energy overes- 
timation. On the other hand, an eventual hybrid event 
of this kind will show a limited, or even below threshold 
FD signal. This leads to contradictory FD and SD en- 
ergy determinations. Moreover, these "golden" events 
would allow identification of i?-hadrons from eventual 
quasi-horizontal neutrino events that are likely to gener- 
ate showers with similar ground signal, but non-negligible 
fluorescence contribution [63]. 

The pertinent question at this point is whether existing 
experiments have already collected events exhibiting the 
characteristics of gluino showers described above. None 
of the ultrahigh energy cosmic ray experiments have thus 
reported such results. It is interesting to note, however, 
that the lO'^ GeV "Centauro" events detected at Mt. 
Chacaltaya [64] might be suggestive of gluino-induced 
showers. In these events, the ratio of hadronic to elec- 



in terms of the Mandclstam variables s = (pa + Pb)'^ — 
{Pc + Pd)"^ and t = ipa - Pcf = {Pb - Pd)'^- Note that 
this squared invariant when viewed from the cm. frame 
reduces to 

I*^ = Vs {Pa X P*) 

= V~s \p:\ \p:\ sine* , (7) 



tromagnetic components is about 50:1, contrary to the 
expectation of dominance of the electromagnetic com- 
ponent in vertical baryon-induced showers. The most 
carefully considered explanation to date is the explosive 
quark matter model [66]. Interestingly, though heavy 
high energy gluinos could also produce such an inverted 
hadronic electromagnetic ratio. This is because the mul- 
tiple low-inelasticity collisions would result in hadronic 
superimposed showers. At detector level 5200 m), 
the electromagnetic component of the sub-showers would 
be mostly filtered out, while the superposed hadronic 
showers would survive. This is because the "low" en- 
ergy (~ 100 TeV) electromagnetic sub-showers induced 
by high energy i?-hadrons would develop faster (being 
quickly quenched by atmospheric losses) than the high 
energy (~ 10^ GeV) electromagnetic subshowers induced 
by ultrahigh energy i?-hadrons. It is also interesting to 
note that this explanation of the Centaruo events does 
not predict any phenomenon one might observe at a col- 
lider experiment, consistent with the null-results from 
UAl [66], UA5 [67], and CDF [68]. If in fact, gluinos are 
guilty of producing the Centauro events, it would consti- 
tute the first evidence of a finely-tuned universe from a 
cosmic ray observation. 
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Appendix 

Consider the process in which two particles of 4- 
momenta Pa and pb and masses Mji and jtitv scatter two 
particles of momenta pc and pd and masses Mji and Mx , 
respectively. Using the total 4-momentum P we define 
the vector 

4 = ec./3M. P'i Pc (5) 
and write the Lorentz-invariant form 



(6) 

I 

where 6* is the scattered angle. Consequently, the for- 
ward direction is defined through the condition lal" = 0. 
In the large s limit where 

s t (2M^ + m% + M]^ - s - t) !v -s'^ t, (8) 

the minimum momentum transfered can be easily ob- 
tained by setting Eq. (6) = and solving for tmin- All in 



I 

1^1'' = St {2M% + m% + M% - s-t) - t {M% - ml) (M^ - M%) - M% (M^ - m%f 
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all, 



Now, combining Eqs. (9) and (10) we obtain 



^min — 



{Ml-Ml){ml-Ml) + s^ 



Ml 



(9) 

ml)/{2^) and 



In the cm. frame, E* = (s 4 

E* = (s + M2 -M|.)/(2yi) [2]. Therefore, the invariant 
quantity {Ec — Ea)/Ea that describes the inelasticity of 
the process reads, 



-^inel 



\t ■ 
Mr 



(11) 



(5 



Ml 



Ml)-{s + Ml) 



is- 



M2) 



Mi 



(10) 



The QCD cross section falls off very rapidly and gets 
negligible for t > Aqcd- Thus, taking Aqcd ~ 1 GeV 
Eq. (11) leads to Ki^ei {Mr/GbY)-^. 
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